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Ageing in systems without detailed balance is studied in bosonic contact and pair-contact 
processes with Levy diffusion. In the ageing regime, the dynamical scaling of the two-time 
correlation function and two-time response function is found and analysed. Exact results for 
non-equilibrium exponents and scaling functions are derived. The behaviour of the fluctuation- 
dissipation ratio is analysed. A passage time from the quasi-stationary regime to the ageing 
regime is defined, in qualitative agreement with kinetic spherical models and p-spin spherical 
glasses. 
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1 Introduction 



One of the paradigmatic example of non-equilibrium critical behaviour is furnished by ageing systems. 
A common way to realise physical ageing is to move a system rapidly out of equilibrium by 'quenching' 
some thermodynamic parameter into a coexistence region of its phase diagramme such that there 
exist several competing thermodynamical states. Alternatively, one may quench a system exactly 
onto a critical point of its stationary state, which is the case we shall study in this paper. From a 
phenomenological point of view, physical ageing may by characterised by the three properties of (i) 
slow, i.e. non-exponential dynamics, (ii) breaking of time-translation-invariance and (iii) dynamical 
scaling. 

Although physical ageing was first recognised and studied in glassy systems, starting from Struik's 
classical experiments [HI H5] , the analysis of many aspects of ageing is more simple in non-glassy 
systems without disorder or Frustrations! Typical examples of this kind are simple ferromagnets, the 
dynamics of which can be characterised in terms of a single time-dependent length scale Lit) ~ t 1 ^, 
which defines the dynamic exponent, see [TTJ [12] for reviews. Their ageing behaviour is conve- 
niently studied through the two-time correlation and response functions of the time-dependent order- 
parameter 4>(t), for which ones expects the scaling behaviour 



C(t,s) = (<f>(t)<f>(s)) - (cf>(t))(cf>(s)) = s- b f c (t/s) (1.1) 

= s-^fnit/s) (1.2) 



R(t, S) =r ^ 



5h{s) 



h=0 



where h(s) is an external magnetic field conjugate to the order-parameter, and t, s are the observation 
and waiting times, respectively. This scaling behaviour should be valid in the ageing regime t,s^$> 
Tmicro an d t — s ^> r micro where T micro is a microscopic reference time. From the asymptotic behaviour 
fc,Ft(y) ~ y~ x c,n/ z f the scaling functions as y — > oo one may define the autocorrelation and 
autoresponse exponents Ac and Xr, respectively. The exponents a, b are known as ageing exponents 
and can be expressed in terms of z and static critical exponents. For reviews, see e.g. [HI [HI |2T| 
251. 



In this paper, we are interested in systems where the underlying dynamics does not satisfy 
detailed balance such that the stationary state is no longer an equilibrium state. Recently, it was 
shown that the same kind of non-equilibrium scaling behaviour (11. 1111.21) as described above for 
ageing ferromagnets applies to this kind of systems, as exemplified for the critical contact process 
[T9| W\\ [8] and the critical non-equilibrium kinetic Ising model [36] . However, in contrast to the non- 
equilibrium critical dynamics of ferromagnets, the ageing exponents a and b need no longer to be 
equal, which for the contact process can be understood as a consequence of its model-specific rapidity- 
reversal-symmetry [8]. Furthermore, there are exactly solvable models such as the so-called 'bosonic' 
contact and pair-contact processes^ [2H1 E3 EE] for which the ageing behaviour could be analysed in 
full detail and the non-equilibrium exponents and scaling functions be calculated exactly [5]. The 
critical behaviour of the 'bosonic' models is different from the more usually studied 'fermionic' ones 
since, although the average total number of particles is conserved at criticality, for long times the 
particles may condense onto a single site which mimics the inhomogeneous growth of bacteria colonies 
and had originally led to the formulation of the BCPD [16"1 |2"%] and BPCPD. All these models show 



2 A careful comparison of the ageing properties of glassy and non-glassy systems and a detailed appreciation of their 
differences was given by Chamon, Cugliandolo and Yoshino [T5] . 

3 The terminology 'bosonic' refers to the property of these models (referred to as BCPD and BPCPD, respectively) 
that an arbitrary number of particles per site is allowed, in contrast to the usual contact or pair-contact processes 
which obey the 'fermionic' constraint of at most one particle per site. 
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ageing only when brought to the critical point (s) of the stationary states, see [23] for a recent review. 
They also have in common that transport of single particles is only through the 'diffusive' hopping to 
nearest-neighbour sites. We wish to investigate the consequences when that local, nearest neighbour 
transport of single particles is replaced by long-range Levy flights, leading to superdiffusive behaviour. 
There are several motivations for such an undertaking: 

1. Levy flights [331 [35] generalise the usually considered local (brownian) random walks in that 
they lead to a generalised central limit theorem for the sum of a large number of independent 
random variables. Heuristically, it appears reasonable to consider the 'particles' which make 
up the basic dynamics of non-equilibrium systems as coarse-grained variables which arise from 
summing over many more microscopic degrees of freedom. It may hence appear natural that 
the effects of long-range jumps between distant sites should be taken into consideration. This 
might become relevant for the description of traces particles in turbulent flows or the spreading 
of epidemics, see [33] and refs there in. 

2. When studying the ordinary contact process with Levy flights, it has been shown that the 
critical behaviour of the stationary state, as well as the relaxation behaviour, is modified 

[271 EH ESI [21 [26]. 

The models under study are defined as follows. Each site of an infinite rf-dimensional hyper-cubic 
lattice may contain an arbitrary positive number of particles. Particles on the same site can undergo 
the reactions 

mA — > (m + k)A with rate /i 

mA — > (m — £)A with rate A (1-3) 

Furthermore, single particles can hop to another site at distance |r|, with a rate 

d r°° 

Dir) = j2nY J ^ e(iq ' r ~ C " qr) (L4) 

where < rj < 2 is a control parameter and c is a non- universal dimensionful constant. The case 
m = 1 (and with k — £ — 1) is called the bosonic contact process with Levy flight (BCPL) and the 
case m = 2 is called the bosonic pair-contact process with Levy flight (BPCPL)o 

In section 2, we write down the closed set of equations of motion for the correlation and response 
functions. In section 3, we discuss the phase diagram, exponents and scaling functions for both 
models. In section 4, we verify that the exact scaling functions found for the bosonic contact process 
can be understood from the theory of local scale-invariance. A discussion of local scale-invariance 
for bosonic pair-contact process with Levy diffusion will be presented in a sequel paper. In section 5, 
we study the passage from the quasi-stationary regime to the ageing regime and identify the passage 
exponent ( which characterises the relevant time-scale t p (s) ~ s^ [47], for the first time in a critical 
system. In section 6, we present our conclusions. Finally, we derive in appendix A a relation for the 
physical interpretation of the two-time correlators, list in appendix B the scaling functions for the 
space-time correlators and responses and present in appendix C the analysis of the passage time in 
the critical spherical model. 

4 The results previously found for the BCPD and BPCPD [28, 37, 5 will be recovered in the limit case 77 — > 2. 
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2 Equations of motion 



Following [37], the master equation is written in the quantum hamiltonian/Liouvillian formulation 
(for reviews, see [42, 26J) as d t \P(t)) = —H\P{t)) where P(t) is the time-dependent state vector and 
the hamiltonian can be expressed in terms of annihilation and creation operators a(x) and a^(x). 
We also define the particle number operator as n(x) = a^(x)a(x). The Hamiltonian reads 



tf = -]T ]T /)(!■) [a(x) a^x + r) - n(x)] 

r^O x 

m 

- A ^[at(x)^a(x) m - JJ(n(x) - i + 1)] (2 . 1} 

X 1 = 1 

m 

- ^[at(x) m+fc a(x) m - JJ(n(x) - i + 1)] - £ Mx)a f (x) 

x i=l x 

For the computation of the response function, we also added an external field which describes the 
spontaneous creation of a single particle with a site-dependent rate h = h(x) on the site x. 

Single-time observables g(t, x) can be obtained from the time-independent quantities by switching 
to the Heisenberg picture. The differential equations for the desired quantities can be obtained by 
using the usual Heisenberg equation of motion d t g = [g,H]. The space-time dependent particle 
density p(£,x) = (at(t, x)a(£,x)) = (a(t,x)) satisfies 



|-<a(*,x)> =^D(r)[(a(t,x + r)) - (a(t,x))] 

9t ^ (2.2) 



- A£(a(t, x) m ) + fjik(a(t, x) m ) + x) 

For the bosonic contact process BCPL, this equation closes for arbitrary values of the rates. However, 
for the bosonic pair-contact process BPCPL, we only find a closed set of equations along the critical 
line defined as 

a = » k - Xe lQ (2.3) 
As we shall see later, along the critical line the space-integrated particle density 

p := J dx p(t, x) = Jdx (a(t, x)) (2.4) 

is conserved, although the microscopic process can change the total number of particles in the system. 
All this is completely analogous to what was already known for the BCPD and BPCPD. 

Throughout this paper, we are interested in the correlation and response functions 

C(t;r) = (a(t,x)a(t,x + r)} — pi 
C{t,s;r) = (a(t,x)a(s,x + r)} — pi (2.5) 
5{a{t,x + r)) 



R{t,s;r) 



h=0 



8h(s, x) 

In these notations, we already anticipate spatial translation-invariance. The single-time correlator is 
determined from the following equations of motion 



-|(a(t, x)a(t, x)> = 2 ^ D(r)[(a(t, x)a(t, x + r)> - (a(t, x) 2 )] 



dt 

r^O 



+ X£[(l + £-2m) (a(t, x) m ) - 2{a(t, x) m+i )] 
-pk[{l -k- 2m)(a(t,x) m ) - 2{a{t, x) m+1 )] 



(2-6) 
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and for x^y 

d 

— (a(t,x)a(t,y)) = ^£>(r)[(a(t,x)a(t,y + r)) + (a(t, x + r)a(t, y)) - (a(, tx)a{, ty))} 

r^O 

- A£[<a(t, x)o(t, y) m ) - (o(t, x) m a(t, y)>] (2.7) 

+ ^[(a(t l x)o(i ! yn - (a(t,x) m a(t,y))] 

In particular, these equations close along the critical line f !2.3j) . The equation of motion of the 
two-time correlator reads 

d 

— (a(t,x)a(s,y)) = V-D(n)[(a(t,x + n)a(s,y)) - (a(t,x)a(s,y))] 

^ ^ (2.8) 

- A£(a(t, x) m a(s, y)) + fik(a{t, x) m a(s, y)) 

together with the initial condition \im t -> s {a(t, x)a(s, y)) = (a(t,x)a(t,y)). 

The relationship of these correlators with the average density and its variance are given for 
diffusive transport by (n(t,x)) = (a(t,x)) and (n(t,x) 2 ) = (a(t,x) 2 ) + (a(t,x)) [37] and more 
generally by |4J, using also the definitions (12.51) 

(n(t, x)n(s, x + r)) - p 2 = C(t, s; r) + p R(t - s, r) (2.9) 

In appendix A, we generalise the proof of these relations to the case at hand. We shall show there 
that the second term in (12.91) is merely generating corrections to the leading scaling behaviour so 
that for our purposes, C(t, s; r) can be interpreted as a connected density- density correlator. 

We begin the analysis of the equations of motion with the particle-density. For the BCPL, as 
well as for the BPCPL with a = 0, one can Fourier-transform eq. (I2.2p . with the result 

J^(f,q)> = -\H ^j^ 1 " ^ + \ <flM) +h(t,ci) (2.10) 

n^O 

V v ' 

= :uj(q) 

which defines the dispersion relation u(q). If we can perform a continuum limit 

w (q) = i-e-cNI" ~ c \\q\\ r > (2.11) 

q^O 

from which we read off the expected dynamical exponent z = rj. 

Similarly, we calculate the response function by applying its definition (12.51) to the equation of 
motion (12. 2\i . This becomes in Fourier space (again, one must set a = for the BPCPL) 

-R{t, s; q) = --w(q)% s; q) + ?-R{t, s; q) + 6{t - s) (2.12) 
at I I 

Consequently, the real-space response function reads 

R(t, s; r) = — i- / d d q e iq r e -X«0(*-«) exp ( a ^~ S A Q(t - s) (2.13) 
(27r) y B V 2 / 



where B is the Brillouin zone and the B-function expresses causality. In particular, for long times 
the autoresponse becomes for t > s 

R(t, s; 0) ~ e§ <*-> ( t -^-) -L f°° d d q e'^ = ^Ii^L§(*-) ( c*—^-) (2.14) 



(27T) d J ^ 7l d / 2 7]T(d/2) 

For the correlation function, we shall assume that spatial translation-invariance holds and use the 
non-connected correlator F(t, s; r) := (a(t, x)a(s, x + r)) and also introduce the control parameter 

a = ^2-b— (2 - 15) 

As initial conditions, we shall use throughout the Poissonian distribution F(0, 0; r) = p§. Generalising 
slightly the calculations performed earlier for the BCPD and BPCPD [371 [5], we find for the BCPL 
the connected correlation function (on the critical line a = 0) 

C(t, s; r) = ap dr &(— !p - r > r ) - a Po f Q dT f exp f _c IMI' ? ~ T ) ~ ir ' 9 ) 

(2.16) 

where the long-time limit t, s S> tmicro 

and t S 3> tmicro 

was also taken. See appendix B for 
the computational details. For the bosonic pair-contact process BPCPL, we first consider the non- 
connected single-time correlation function F(t, r) which satisfies the equation 



F(t,r) =p 2 + a [ drF(r,0) / 

JO J B 



t r Ad r _ 

w(q)(t--T) iq-r 



d d q 



b (2vr) d 



/•< (2 ' 17) 
= p 2 + a / dr F(t, 0) 6(t - r, r) 
Jo 

For r = 0, this is a Volterra integral equation for F(t, 0) which may be solved by Laplace transfor- 
mation as in ([37], [5]). The F(t,r) is then directly obtained from (12.171) . Having found this, the 
two-time correlator is given by 

F(t, s,r)=p 2 + a^ drF(r, 0)6 - r, rj (2.18) 

The results of eqs. ( 12. 16112. 17)2.180 will form the basis of the subsequent analysis. The essential 
difference with respect to the BCPD and BPCPD is the different form of the dispersion u)(q). 



3 Exact solution 
3.1 Phase diagramme 

Since the total number of particles is conserved on average, information on the critical behaviour 
comes by analysing the variance (n(t,r) 2 ). The resulting phase diagramme is presented in figure 
CD There are two distinct phases separated by the critical line eq. ( 12.30 . namely (i) an absorbing 
phase for A > ji with a vanishing particle-density in the steady-state, and (ii) an active phase for 
A < /i, where the particle-density diverges for large times. Along the critical line, the total mean 
particle-density is constant. The behaviour of the model along the critical line, described by varying 
the control parameter a eq. (12.150 is as follows: 
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A 




(2) 



W (b) 



Figure 1: Schematic phase-diagrammes for d 7^ of (a) the BCPL and the BPCPL in d < rj and (b) 
the BPCPL in d > rj. The absorbing region 1 for A > /i, where the particle-density vanishes at large 
time, is separated by the critical line eq. ( I2.3P from an active region for A < /z, where the particle- 
density diverges at large times. Along the critical line, the average time- dependent particle-density 
remains constant. 



• For the bosonic contact process BCPL, one has the same critical behaviour for all values of 
a. On the other hand, the value of the space dimension d is very important. If d < rj, one 
has clustering: for long times, the particles are on average redistributed such as only a few, 
spontaneously selected, lattice sites contain particles while the others become essentially empty. 
For d > r), the long-time particle-density is spatially homogeneous. 

We can interpreted this result by using a generalised Polya theorem [39] , which states that for 
d < rj the Levy random walk is recurrent and hence cannot homogenise the system. On the 
other hand, when d > rj, there is a non-zero probability that a Levy random walk does not 
return to its starting point. This is enough to make the system spatially homogeneous. 

• For the bosonic pair-contact process BPCPL, when d < rj, there is always clustering. On 
the other hand, for d > rj, there is a multi-critical point at a = etc such that the system is 
homogeneous for a < etc and that one has clustering for a > etc- The critical point of this 
clustering transition is given by 



From this expression, one sees that ace — > oo if d — > rj [2lT [6]. 
3.2 Ageing exponents and scaling functions 

Using the explicit expressions from section 2, we can perform the long-time limit and by comparing 
with the scaling forms eqs. (11. Ill . 21) we can proceed to extract the non-equilibrium critical exponents, 
which are listed in table [1] and then the scaling functions which, up to a normalisation factor, are 
given in table EJ For the details of the calculation, see appendix B. 

Some comments are in order. 

• If d > rj and if furthermore a < ac in the BPCPL, the anticipated dynamical scaling (ll.l|1.2l) 
holds true for both the BCPL and the BPCPL. Furthermore, since the exponents and the form 




(3.1) 
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Bosonic 
contact process 


Bosonic pair-contact process 


a < otc 


a = otc 


a 


d_i 

V 




n 


b 


a - 1 

V 


V 


if r] < d <2r] 

1-2 if d > 2ri 
n 1 


Xr 


d 


d 


d 


Ac 


d 


d 


d 


z 


7] 


V 


V 



Table 1: Ageing exponents of the critical BCPL and BPCPL in different regimes. The results for 
the BCPL hold true for any dimension d, while for the BPCPL they only apply if d > rj. 





fn(y) 


fc(y) 


Bosonic contact process 


(y-i)- d/v 


(y-l)-d/V+l _ (y + l)-<i/V+l 


Bosonic 
pair contact 
process 


a < ac 


d > r] 




( y _ l)-«t/«7+l _ (y + l)-d/ V +l 


a = ac 


rj < d < 2r) 


(y-l)-*/l 




d>2r] 


(y-l)-<l/v 


2(d-2rj)(y- l) 1 ^ 
-r]{y - l) 2 - d A? + rj(y + lf- d h 



Table 2: Scaling functions of the autoresponse and autocorrelation function of the critical bosonic 
contact and pair-contact processes. They are given up to a multiplicative factor. 

of the scaling functions agree, the two models are in the same universality classE] 

• For d < rj, we still find a dynamical scaling behaviour in the BCPL. We observe that the 
exponents a = b now become negative. This reflects the eventual condensation of the particles 
which means that the fluctuations diverge for large times. 

For the BPCPL however, for d < rj or more generally a > ac, no dynamical scaling is found. 
Rather, the variance diverges exponentially with time. The notion of ageing as defined in the 
introduction no longer applies here. 

• Finally, for the BPCPL with d > rj and at the tricritical point a = ac, the dynamical scaling 
f ll.l|ll.2l) holds true. The values of the non-equilibrium exponents are distinct from those of the 
universality class discussed above and in particular, we note that a ^ b although still Ac = Xr- 
In this respect, the absence of detailed balance has led in this case to an ageing behaviour 
intrinsically different from the non-equilibrium dynamics of a ferromagnet (where a = b as well 
as Xc = Xr always hold true). 

• The response function R = R(t — s; r) has a particularly simple form which actually obeys 
time-translation-invariance. However, since time-translation-invariance is not satisfied for the 
two-time correlators, the notion of ageing is still applicable. 

If we identify the upper critical dimension as d* = rj, we see that our results can be mapped onto 
those of the BCPD and BPCPD if one replaces d/rj = d/d* f-> d/2. 

5 In many respects, the behaviour of the BCPL found here is analogous to the ageing behaviour of a simple 
ferromagnet quenched onto its critical point. 
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4 Local scale- invar iance 



We now inquire whether the forms of the scaling functions derived in the previous section can be 
understood from some larger dynamical symmetry than mere dynamical scaling. 

Indeed, for simple ferromagnets undergoing phase-ordering kinetics after a quench from a totally 
disordered initial state to the ordered phase with a temperature T < T c , an extension of dynamical 
scaling [H] to a local group of time-dependent scale-transformation related to a subgroup of the 
Schrodinger group (without time translations) has been found [231 125] . In particular, it has been 
understood how to use the necessarily projective representations of non-semisimple groups as the 
Schrodinger group in order to analyse the dynamical symmetries of the stochastic Langevin equations 
underlying these phenomena [40] . The essential tool in this kind of analysis are the celebrated 
Bargmann superselection rules [3]. Then both response as well as correlation functions can be 
explicitly calculated and the results have been successfully tested in numerous models, see |25] for a 
review. Since the Schrodinger group and its subgroups only apply to physical situations where the 
dynamical exponent z = 2, a generalisation to arbitrary values of z must be sought. 



4.1 Background 



In this section, we shall use the BCPL with z = rj < 2 as an analytically treatable test case for such 
a possible extension. In order to do so, we shall first reformulate the problem as a non-equilibrium 
field-theory using the Janssen-de Dominicis theory pj2, [30J, [29] . Starting from the master equation 
with the reaction rates eqs. ( 11. 3)1. 4ft . the creation and annihilation operators become related in the 
continuum limit to the order-parameter field and a conjugate response field 



<f)(t,r) :=a(t,r)-p 
0(i,r) :=at(t,r)-l 

such that (0(t,r)) = (<f>(t,r)) = 0. The action associated to the critical BCPL reads 



J[<f>A] 



dR du 



0(2A^ n -A^ 2 )0-/i0 2 (0 + po) 



Jol 



(4.1) 
(4.2) 



(4.3) 



where we have suppressed the arguments of 4>{u, R) and <p(u, R) under the integrals. The dimen- 
sionful 'mass' M. is related to a generalised diffusion constant. The action is decomposed into two 
parts: a so-called 'deterministic' (noiseless) part 



Jo 



dR / du 



<f)(2Md u - A"/ 2 )0 



and a 'noise' part 



dR / du 



/Z0 2 (0 + p O ) 



(4.4) 



(4.5) 



Here the non-integral power of the Laplacian has to be interpreted as a fractional derivative [H [9]. 
In this formalism, correlators and response functions are found as follows: 



tni fl j ■ ■ ■ Tn) 

R(t,s;r,r') 



(0i(ti,ri) . . .(f) n (t n ,r n )) 
(0(t,r)0( S ,r')) 



(4.6) 
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where the average of any observable A is given by the functional integral (A) = f T>(f>D(f)A[(f), <f>]e . 
As we have shown earlier for the case z = 2 [ID], one considers first the 'deterministic' part of the 
action, obtained from the variation of j7o with respect to <fi. The associated equation of motion 

S<f) = (2Md t - A r ' /2 ) = (4.7) 

has a dynamical symmetry under the infinitesimal generators of local scale-transformations [U [9] 

X_! := —d t time-translation 

1 x 
X Q := —tdt (r- d r ) dilatation 

77 77 



Xi:= .e dt .^±H t .^l-,.l t{r . dr 



rj rf " rj 



generalised Schrodinger 



—27(2 — 77) (r • 9 r )V r r? — 7(2 —rj)(d — ^)V r r ' transformation 

Y_h := — d n space rotation 

Y® /v+1 := -td n - ^ 1 —r i V 2 - r > - 777(2 - v)d r ^r V generalised Galilei (4.8) 

transformation 
R(^) ■= -rM = n d r . - Tjd ri rotation 

with i,j — 1, . . . , d, 7 is a dimensionful constant and x,£ are scaling dimensions. The fractional 
derivatives V° with a6l are defined in [U,[9]. From the commutation relations [HE] 



[A n ,X n /] — (n — n )X n+n i 



n+m 

[Y$>,BM] = y® (4.9) 

[X nt RfM]=0 ; [Yi k \R^}=0 i£k?i,j 

it is clear that the explicitly specified generators span the complete algebraic structure. Furthermore 
[S, X] — with all generators X of the above list (14. 8p . with only two exceptions, namely 

[S,X ] = -s 

[S,X 1 ] = -2tS~?Y (2(x + 0-( V -2 + d)-^Lrj 2 (2-ri)) (4.10) 

Hence the infinitesimal generators ( 14.81) map a solution of the 'deterministic' equation S(f) = onto 
another solution, provided <j) the scaling dimensions x, £ of <ft satisfy the condition 

l + {=2z 2 ±f? + ^ (2 _ )|) (4n) 

In order to be able to analyse the full theory, we first define the so-called 'deterministic averages' 
(A) := jVcffDcf) A[cf), 0]e~" 7 °^''^ [40]. Consideration of the iterated commutators of the generators 
Ym then shows that these deterministic averages obey generalised Bargmann superselection rules 

[HE] 

0, unless n = m (4-12) 





Now all building blocks are provided that we can adapt the treatment given earlier for z = 2 [40] to 
the present case. 
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4.2 Application to the BCPL 

1. In order to compute a response function, one treats the 'noise' part of the action as a 
perturbation. Using the relations (14. 6[) and (14.121) . we have 

R{t, s; r, r') = (^>(t, r)0(s, Y ') e -M4>M^ 

= (0(t,r)0( S ,r')) o 

= : ^, S )^^)(^-^) (4.13) 

and where we used the explicit form as determined from the covariance of the 'noiseless' response 
function. This shows that the response function does not depend explicitly on the noise and may 
be found directly from the symmetries (14.81) of the 'deterministic' part alone. In particular, the 
form of the autoresponse function R(t, s) = R(t, s; r, r) is determined by its covariance under 
the generator X\ and reads 

/ , \ l+a'-\ R /z s, \ -l-o' 

-m-i f--i 



R(t, s) = s-"" 1 ( - 1 (--11 (4.14) 

This result is in agreement with the scaling form (table [2]) derived in section 3 with the expo- 
nents a' = a = d/i] — 1, X R = d. 

2. We have for the space-time correlator, generalising the treatment of [6] to z = rj < 2, and 
following p[9] 

C(t, s; r, r') = (^(t, r)0(s, r ')e~ Jb[ ^ ] } (4.15) 

Expanding exponentials and using the Bargmann superselection rule, we find that the correlator 
is the sum of two terms Ci(t, s; r, r') and C 2 (t, s; r, r') 

Ci(t, s; r, r') = -2fip J dR J du ($(t, r)0(s, r')<p 2 (u, R)^ (4.16) 

and 

C 2 {t, s; r, r') = J dRdR' J dudu'{<f>{t, r)(f){s, r')T{u, R)T{u', Rf)) (4.17) 

where T := 2 is a composite field with a scaling dimension xx = 2x + x. Furthermore, both 
terms of the correlator can be factorised as a product of response function jl] such that 

C 1 (t,s;r,r') = -2//p J dR J (0(t, r)${u, R)) ■ (<f>{s, r')4>{u, R)) 

= -2/ipos- 1 - 2 " / du / d d Ru- 3 - 3a+XR/ri ( — -l] 

Jo Ud \su J \u ) 

x (L\ ^ A) ( (r-R) \ {a2A) ( (r>-R) \ 

\su) \sVi{\-uyh) \s^(t/s-u)^J { } 
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and 



C 2 (t, s; r, r') = 2/x 2 j dR J du (0(t, r)T(u, R)) Q ■ (<f>{s, r')T(u, R)) 
= s 1+ ^^2/i 2 / du / d d i?M^^ 1 



JO JR d V' 5 " 

"is- 1 L) < 4 - 19 > 



X 



Using a = — 1, a dimensional analysis shows that 

C n {t, s; r, r') = s^^F^t/s, rs" 1 ^, r's' 1 '^) (4.20) 

with n — 1,2 and _F n are scaling functions. It follows that in the scaling regime s — > oo 
(with the other scaling variables fixed) the term C% is negligible. It remains to see that C% is 
compatible with the expressions given in Table [2j Using the identity [I], 

/ dRdR' (cR + a) T^fr) ( dR + 6 ) g ^ 

= I T^Yd I ^/ ak+ ' lbq 9(ck,dq)\k\^^ (4.21) 

To recover the expression given in eq (14. 181) . we require g(R, R') = 5{R — R') and its Fourier 
transform reads (27r) d 5(k + q) . Then, we can easily find, with fli — ft% — f3 and a.\ = a%, 

d(t,s;r,r') = -2/ipos f du u- 2 ~ 2a+2XR/v (t/s - u) -1 ~ a+Aj?/r? (l - u)- 1 -^/" 

Jo 

x f d k lfe|2/3 e i(r-T-0-fe-2a 1 |fc["(t+ S -2n S ) ( 4 22 ) 

Since a = d/rj — 1 and A# = d, the full space-time correlator reads 

s; r, r') = -2 m f d^ / V^W^^e*^* (4.23) 

Jo Jr* ( 27r ) 

where we identify (3 = 0, aj = ^ = f > and M. = — to match to the expression (I2.16p . 

Therefore, although the field-theory of BCPL is not free, the structure of its 'deterministic' part 
eq ( 14.7ft is simple enough to explain the exact results as a manifestation of the local- scale-invariance. 



5 Passage time towards the ageing regime 

In the previous section, we have analysed the scaling behaviour of the BCPL and BPCPL. Now, 
we analyse how from an initial state this scaling regime may be reached. Since the BCPL and the 
BPCPL with a < ac are in the same universality class, we shall limit ourselves in what follows to 
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Figure 2: Autocorrelator C(t, s) of the critical 3D BCPD (77 = 2) as a function of t — s, for several 
values of s. Full lines: complete autocorrelator. Dash-dotted line: stationary correlator eq. ( 15. 2ft . 
Symbols: scaling contribution C ag emg(t,s) = s~ b fcit/ s). 



the BCPL or even to the BCPD. Our analysis generalises the previous study of Zippold et al. |47j . 
who examined this question for the low-temperature spherical model and the spherical spin glass, in 
that we look at the critical case. 

For a qualitative overview, we show in figure [21 the autocorrelation function C(t, s) as a function 
of time separation t — s. For small time-separations, the autocorrelator remains close to the quasi- 
stationary, time-translationally-invariant form implied by the Poissonian initial conditions, before 
crossing over to the scaling regime with ageing behaviour when t — s becomes large. In order to be 
able to define a precise measure of this cross-over, we consider the following relative error 



6C{t, s) 



C(t,s)-C stat (t,s) 



(5.1) 



C a geing(tj Cstat(^) 

where C stat (t, s) is the quasi- stationary correlation function 

/»00 /*7T 

C sta t{t, s) = ap dr / dqq d ~ l exp (w(g) r) (5.2) 

Jo 

and C a geing(t, s) = s~ b fc(t/s) is the autocorrelator in the ageing regime, see eq. (12.161) . Since < 
SC(t, s) < 1, we arbitrarily consider the system to be in the stationary regime when 5C(t, s) < 10% 
and to be in the ageing regime when 5C(t,s) > 90%. In figure El this relative error is plotted as 
a function of the separation time for different values of s. We observe that the passage from the 
quasi- stationary regime to the ageing regime occurs at relatively well-defined transition-times. 

Two transition times will be defined: 

1. the system leaves the quasi-stationary regime at the time scale r stat (s) which is defined by 

SC(T stat (s) + s,s) = 0.1 
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d=3, T]=2 d=5, ri=2 




Figure 3: Relative error 8C(t,s) as a function of time separation t — s for different values of the 
waiting time s in the BCPD. Left panel: d = 3, right panel: d = 5. 



2. the system enters the ageing regime at the time scale r ag eing(s) which is defined by 

<5C(Yageing(,s) + S, s) = 0.9 

These transition times are shown as a function of the waiting time s for several dimensions in figure HI 
We find for sufficiently large s the asymptotic behaviour 

1~ageing{t > s ) = As (5-3) 

where A depends on dimension and the passage exponent £ < 1. In principle, ( should depend on 
both the dimension d and the Levy parameter 77. However, as illustrated in figure [5j £ apparently 
depends merely on the ratio d/r/. Indeed, we obtain a neat collapse of the entire curves r a g ein g(s) for 
different values of d and 77, but with the same value of d/i]. Therefore, it is enough to look at models 
with a fixed ratio of d/rj and for that reason most of our calculations were performed for the BCPD. 



The dependence of ( on the dimension d is shown for rj = 2 in figure [5] and in table [3) Analogously 
to the study of Zippold et al. [47] for spherical ferromagnets and spherical spin glasses quenched to 
T < T c , our results for the critical passage exponent indicate that ( < 1. In table [3] we list some 
known values of ( (the values for the critical spherical ferromagnet are obtained in appendix C). 
In conclusion, we find that the same qualitative behaviour for the passage time between the quasi- 
stationary and the ageing regime, r a g eing (s) ~ with < ( < 1 apparently holds true both for 
critical as well as for non-critical systems. This form of the passage time is one of the ingredients 
for the derivation of admissible scaling forms of the two-time autocorrelator, which gives C(t, s) = 
C st (t - s) + C age (t, s), where [TJ 



Cage(^) s) — C 



h{s) 



h{t) = h exp 



1 t 1 -^ - 1 



A 1 



V J 



(5.4) 
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Figure 4: Transition times r stat (s) (lower curves) and r ageill gi 
s of the BCPD, in 3D (left panel) and in 5D (right panel). 
7.7 



s) (upper curves) as a function of 
The full lines give the linear fits: 



7~agcing \ & 



s o.3i f or ^ = 3 anc i r ageing (s) = 


12.6 • s - 53 


for d = 


5. 




model 


condition 


d 


c 


Ref. 


XY 


T = 


1 


2/3 


[22] 


spherical 


T <T C 


> 2 


4/(d + 2) 


02] 


p-spin spherical glass p = 2 


T<T C 




4/5 




p = 3 


T<T C 




« 0.68 


[32J 


BCPD 


\ = fi 


3 


0.31(1) 








3.5 


0.40(1) 








4 


0.47(1) 








5 


0.53(1) 








6 


0.58(1) 




spherical 


T = T C 


3 


0.33(1) 








5 


0.40(1) 





Table 3: Values of the passage exponent £ defined in eq. (15.31) for several models, either in the 
coexistence phase or at criticality. 



where C is a scaling function, [l is a free parameter related to ( and ho and A are normalisation 
constants. However, while for quenches for T < T c the available evidence suggests that ( should 
decrease with d, see table [3] [47] , we find the opposite tendency for critical quenches, see figure [5] and 
table [3j 

The passage towards the ageing regime can still be understood heuristically in a different way. 
At equilibrium, the fluctuation-dissipation theorem (FDT) describes the expected size of the fluctu- 
ations, given the response to an external perturbation. If the FDT is broken in ageing phenomena, 
the inverse fluctuation-dissipation ratio (FDR) describes by how much this expectation is larger than 
the actually found result, for a given value of the temperature T [33]. Since for the BCPL, the ageing 
exponents a = b are equal, we can define an analogue of a fluctuation-dissipation ratio 



X(t,s) 



dC { 

ds 



s,s) R(t,s) 



R(s,s) 



(5.5) 



15 




s 

Figure 5: Dependence of the transition time r age i ng (,s) on the ratio d/r]. The upper set of curves 
corresponds to d/rj = 4 and the lower one corresponds to d/rj = 3. 




Figure 6: Exponent ( against dimension d for the BCPD (where r\ = 2 was used). If d is replaced 
by rjd/2, the results are also valid of the BCPL. 



but where now the initial FDR plays the role the temperature in system with an equilibrium sta- 
tionary state. From our exact solution of the BCPL, we find, in the ageing regime 



-^ageing j Sj 



1 



(5.6) 



In figure [7J we show X(t, s) and ^ ageing (t, s) as a function of the time separation t — s for different 
waiting times s. We can distinguish three different regimes: 

1. a quasi-stationary regime with microscopic relaxation for t — s r stat (s) 

2. a non-analytic transition regime for r stat (s) <C t — s r ageing (s) 

3. the ageing regime when t — s ^> T agoing (s). 
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The bold circles in figure [7] indicate the transitions between these regimes. We see that the FDR is 
still very close to unity when the quasi-stationary regime is left and not too far from its limit value 
Xoo = \ upon entering the ageing regime. 



d=3, r|=2 




t-S 

Figure 7: Fluctuation-dissipation ratio X(t, s) (full lines) and the ageing limit -X" ageing (t, s) (symbols) 
as a function of the time separation t — s for different value of s in the 3D BCPD. 



6 Conclusions 

In analysing the non-equilibrium dynamical scaling of two 'bosonic' particle-reaction models with 
Levy-flight transport of individual particles, characterised by the parameter r], instead of diffusive 
transport, we have obtained the following results: 

1. We have defined the 'bosonic' contact and pair-contact processes BCPL and BPCPL. Because 
of the absence of a 'fermionic' constraint limiting the number of particles per site, the models 
are exactly solvable, either in the full parameter space as the BCPL or along the critical line as 
the BPCPL. Since the global particle number is conserved on the critical line, information on 
the nature of the phase transition comes from analysing the variance of the particle distribution. 
This shows that the BCPL and BPCPL may give rise to a condensation of the particles on 
essentially a single lattice site. 

2. Starting from an initial condition far from the stationary state, the two-time correlation and 
response functions obey the same kind of dynamical scaling and ageing, see eq. (11.1111.21) as 
already found in other systems without detailed balance. The explicit results for the non- 
equilibrium exponents and the scaling functions are given in tables [Hand [2] for arbitrary values 
< 7] < 2. 

For i] = 2 our results include those of the diffusive case (models BCPD and BPCPD [2"g j [57 ] 15]) 
as special cases. 
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3. We have used these explicit results to test a formulation of local scale-invariance for a dy- 
namical exponent z = i] < 2. Extensions of global dynamical scaling to a more local form in 
non-equilibrium systems must be capable of treating the dynamical symmetries of stochastic 
Langevin equations, which is done by attempting to perturb around the noiseless, 'determinis- 
tic' part of the model, which in the case of the BCPL studied here takes the simple Markovian 
form (2Md t - A^ 2 )0 = 0, where M is a dimensionful constant. Local scale-invariance shows 
how generalised Bargmann superselection rules can be derived for arbitrary values of the dy- 
namical exponent such that the perturbation series naturally truncates. Our analysis shows 
how the procedure may be extended to systems without detailed balance and in this sense 
generalises earlier studies performed on long-ranged ferromagnets [7[ [18] . 

4. Finally, we used the exact solution to follow precisely the passage between the quasi-stationary 
regime and the late-time ageing regime. In analogy with earlier studies on this passage in non- 
critical systems (43, [22], our results suggest that the passage time depends on the waiting time s 
also for critical ferromagnets, as well as for systems without detailed balances, as r age i ng (s) ~ s^, 
with < C < 1- 

All in all, the explicit study of the exactly solvable BCPL and BPCPL processes has permitted us to 
control and to confirm several assumptions, usually admitted in studies of non-equilibrium dynamical 
scaling, and hence to extend their range of applicability. 



Appendix A. Proof of eq. (12.91) 



We generalise Baumann's proof pi], valid for diffusive motion of single particles, to the case of 
long-range jumps. We work in the quantum hamiltonian/Liouvillian formalism. By the property 
(s\a\x) = (s\, we can verify that along the critical line 

d t (s\a(x) = (s\[a(x),H] = ^(s|£>(n) ((a(x + n) - a(x)) (Al) 

The Fourier transform is defined as (s|a(k) = e _lk ' a '(s|a(x), therefore 

(s\a(t, k) = exp (~uj(k)(t - s)J (s\a(s, k) (A2) 

where w(k) is the dispersion relation. Using Fourier's theorem (s|a(x) = (2ii)^ d f B dk e~ tk ' x (s\a(\tx) , 
we have for all times s < t 

* h(2V) V 2 ^ (A3) 

= ^R(t - s;x - y)(s\a(s,y) 
y 

where we have used the expression ( 12. 13f) of the response function. Next, we apply ( 1A3I) to the 
two-point correlation function 

F(t, s; r) = (a(t, x)a(s, x + r)) = ^ - s; x - y)F(s, r + x - y) (A4) 
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We also calculate the density-density correlator and we find 

{n(t, x)n(s, x + r)) = (a'(t, x)a(i, x)a^(s, x + r)a(s, x + r)) 

= — s; x — y)(a(t, x)a^(s, x + r)a(s, x + r)) 

^ ( A5 ) 
= 2J R(t - s; x - y)F(s, r + x - y) + R(t - s; r)p 

y 

= F{t,s;r) + R{t-s;r)p 

where we have used the time-dependent commutator of a and a* and also applied (1A4I) . Going back 
to the definition ( 12. 5p the assertion ( 12.91) follows. 

Finally, using the scaling forms (11.1111.21) together with the exponents of table (U we explicitly see 
that the second term in eq. ( 12.91) becomes negligible compared to the first one in all cases. 



Appendix B. Analysis of the correlators 

We discuss those aspects in the calculation of the correlators which go beyond the local BCPD and 
BPCPD models [371 [5]. First, we reconsider the Green's function b(t, r), see eq. ( 12. lip . In the ageing 
regime (t, s ^> 1 et t — s ^> 1), only small values of q will contribute. Therefore, 

From this, we identify the scaling variable £ = rt' 1 ^ and read off the dynamical exponent z = rj. 

For the explicit calculation of the correlator in the BCPL, we merely have to evaluate eq. ( 12. 1611 . 
For the autocorrelator, we find 

• Case d — rj: 



C{t, s; 0) = ap B dr r = ap B In (j^jj = fc(t/s) (B2) 

where B = ^,^1^ and we read off the exponents 6 = and Ac = d. 
Case d 1]: 

nn n, "Po^o ( ( t + s\^ d/ri (l ■ * x 1 

C(t,s;0) 



1 - d/rj 



a Po B o a i-d/ n ( 1 N 



s 1 '^ 1 - (B3) 



and in the second line we looked at the asymptotic form if t/s — » 00. We read off the exponent 
b = —1 + d/r] and A^ = d. 
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The space-time-dependent scaling functions can be found similarly. We merely quote the result: 

t-sV-"" „ I T \/i + S \ 1 -V' N 



where the function C_i reads 

C_ 1 (rr 1/v ) = t d/r >B 1 [ d d q e iq - r e-^ n (clq^y 1 (B5) 

For the BPCPL, using Laplace transform on eq (I2.17P and applying a similar analysis than in [21], we 
evaluate F(t, 0) in the different cases and then use the eq (12 . 1 8[) to compute the connected correlator. 
The most interesting cases are a < ac and d> r\ 



1. a < ac and d > rj: here, F(t, 0) = a °_ a and the scaling function reads 

aacPpBo 

(a - a c ) (f - 2) 



fc(y) = a J p2 ° B ° r [(y + i) 1 -^ - (y - l) 1 "^] (B6) 



2. a = ac'- 

„2 



For 77 < d < 2?7, F(t, 0) = ac , Bo | r(1 ^ /??) | r(d/??) ^ /r? 1 and the scaling function is given by 

fM _ pI ( y + i Y d/ri F (d d d 2 

jew |r(i - d/»7)|r(d/77)d/»7 v 2 J 21 v^'^'^ y + i 



For d > 2i] 

A 2 2- d / r >+ 2 (d-r]) (d-2rj) 



HV) = A ^ lr ,J?"\^ ^ [2(d - 2^) (y - l) 1 ^ - ^ - l) 2 ^ + r?(y + 1) 2 ~^ 



(B8) 

where A 2 = (2 J B ^(2co(q))-^ \ 

3. a > ac or d < rj: F(t, 0) shows an exponential behaviour and we do not have scaling behaviour 
in these cases. 

The space-time-dependent scaling functions can be found from eq (12.181) . 

1. a < ac and d > rj: 

r \ 1 1 — s\ _ / r \ / 1 + s x 



C(t,»;r) = ap„B„ | d ( j ^— J - £7-, | | — | | (B9) 

2. a = ac-: 

For 7] < d < 2r) 



t , 1 \ n 



n=0 
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where C n is defined as 



C n (rt- 1/ ^)=t d/r 'B 1 [ d d q e [qr e'^ 1 (c\q\ n ) n (Bll) 

For d > 2i] 

C^ r) ^B^{^) g {{ LlJr (B12) 
Appendix C. Critical spherical model 

We analyse the passage towards the ageing regime in the critical spherical model. As introduced by 
Berlin and Kac [10J, its spins S x e 1R are real variables on a hypercubic lattice A in d dimensions 
which obey the so-called spherical constraint 

where Af is the total number of spins. The dynamics is given by the stochastic Langevin equation 

d t S(t, r) = V 2 r S(t, r) + i fl \ng(t)) S(t, r) + V (t, r) (C2) 



2 \dt 

where rj(t, r) is a centred gaussian noise such that {rj(t, r)rj(t', r')) = 2T c 5{t — t')5{r — r') and g(t) acts 
as a Lagrange multiplier to enforce the spherical constraint, see below. We only quote from the well- 
known solution [2U] those results we need. The autocorrelation reads (we use an infinite-temperature 
initial state throughout) 

C(t, s) = m r )S (s, r) ) = (/ (*±i) + 2T C f «f - ,) (CS) 



(C4) 



where g(t) is given as the solution of the Volterra integral equation 



g(t) = f(t) + 2T c f &t'f(t-t>)g(t' 
Jo 



and f(t) reads for short-range interactions (Iq is a modified Bessel function) 



B 



Finally, the critical temperature T c is given by 



T c (d)= (2^°°dt(e- 4 %(4t))^ (C6) 



We also recall the scaling functions of the autocorrelator in the ageing regime 

M = T^^^^-l)^ 2<,<4, 
fc(y) = T c ^^((y-iy-^-(y + ir^) d>4 (C7) 
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s 

Figure 8: Passage time towards the ageing regime in the spherical model at T = T c ,as a function 
of the waiting time s, for different values of the dimension d. The full lines give the linear fits: 
Tageing(s) = 4.9 • s 032 for d = 3 and r agcing (s) = 5.2 • s a4 ° for d = 5. 



As in the main text, we define the passage times between the quasi-stationary and the ageing 
regimes. In figure [HI we show r agoing (s) ~ / as a function of s. The full curves indicate the 
numerical solution of the Volterra integral equation flC4j) whereas the symbols were obtained from 
the asymptotic forms flC7j) . Hence r age i ng (s) shows the same kind of qualitative behaviour as the 
BCPL and the resulting values of ( are included in table [31 
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